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ABSTRACT 

We discuss chemical enrichments of ^4000 SDSS early-type galaxies using as 
tracers a large variety of element abundance ratios, namely [C/Fe], [N/Fe], [O/Fc], 
[Mg/Fe], [Ca/Fe l and [Ti/Fej. We utilise the stellar popu lation models of absorption 
line indices from lThomas. Maraston fe Johanssonl (|2011bf) which are based on the MF 
LES stellar library. We confirm previous results of increasing age, [Z/H] and [O/Fe] 
ratios (most often represented by [cc/Fe] in the literature) with velocity dispersion. 
We further derive identical correlations with velocity dispersion for the abundance 
ratios [O/Fe], [Mg/Fe] and [C/Fe], implying that C/Mg and C/O are close to solar 
values. This sets a lower limit on the formation time-scales and star-burst components 
of early-type galaxies to ~0.4 Gyr, which is the lifetime of a 3M© star, since the full 
C enrichment must be reached. [N/Fe] correlates with velocity dispersion, but offset 
to lower values and with a steeper slope compared to the other element ratios. We do 
not find any environmental dependencies for the abundances of C and N, contrary to 
previous reports in the literature. [Fe/H] does not correlate with velocity dispersion 
over the entire parameter range covered, but for fixed age we find a steep trend for the 
[Fe/H] -a relation. This trend is weaker than the analogous for total metallicity (which 
also shows steeper trends at fixed age) owing to the lower Fe contribution from SN 
la for more massive early-type galaxies. We find [Ca/Fe] ratios that are close to solar 
values over the entire velocity dispersion range covered. Tentative, due to large scatter, 
the results for [Ti/Fe] indicate that Ti follows the trends of Ca. This implies a signifi- 
cant contribution from SN la to the enrichment of heavy a-elements and puts strong 
constraints on supernova nucleosynthesis and models of galactic chemical evolution. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
abundances 



1 INTRODUCTION 

The chemical compositions of stellar atmospheres are tra- 
cers of the element abundances of the parent gas clouds 
forming the stars throughout the formation history of a ga- 
laxy. Some elements a re also affected by dredge-up durin g 
stellar evolution (e.g. ISweigart. Greggio fe Renzin H ll989l ). 
Stellar populations are therefore a powerful tool to extract 
information on chemical evolution in the Universe. Ele- 
ment abundances can be directly determined for individu- 
al stars of resolved stellar populations in the Milky Way 
or in nearby dwarf galaxies, using absor ption lines measu- 
red in high resolution stellar spectra (e.g. lEdvardsson et ail 
19931: iFuhrmannll 19981: iBensbv. Feltzing fe Lundstromll2004l : 



Feltzing. Primas fe Johnson! 120091 ; IBensbv et all l201oT) . 
Light averaged spectra must instead be used for determining 
element abundances of distant unresolved stellar populati- 
ons. The absorption features of such spectra are sensitive to 



multiple elements due to velocity dispersion b roadening. The 
Lick system of absorpti on line indices (e.g. IWorthev et al.l 
11994 iTrager et al.lll998r ) have been frequently used for mea- 
suring 25 prominent absorption features in galaxy spectra. 



Elements are produced in stellar nucleosynthesis besi- 
des the primordial nucleosynthesis of H and He. The che- 
mical enrichment of stellar populations depends on the star 
formation history, initial mass function, fraction of explo- 
ding supernovae etc. The chemical pattern of the parent gas 
clouds will be carried on to new stellar generations. Thus 
the chemical enrichment of stellar populations is also affec- 
ted by mechanisms affecting the interstellar medium (ISM) 
such as the efficiency of stellar winds to mix newly syn- 
thesised elements with the ISM, efficiency of outflow from 
galactic winds to remove enriched g as, inflow of less enri- 
ched gas from gas reservoirs etc. (e.g. iMatteucci fe Francois! 
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ll989l ; lMatteuccllll994 ). Chemical enrichment sets stringent 
constraints on galaxy formation and evolution. 

Studies beginning in the late 1970's have revealed 
non-solar abundance rati os for the stella r populations of 
early-type galaxies ( e.g. O'Conncll 1976; Peterson 1 19761 ; 
Burstein et ail Il984; IWorthev, Faber fe Gonzalez! Il992l ; 



Davies. Sadler fc Peletierlll993l ; ISurma fc Benderlll995l T in- 
dicating different chemical enrichment histories. This trig- 
gered more detailed investigations showing that the ra- 
tio between a-elements and Fe-peak elements increases 
with increas i ng galaxy mass for early-type galaxies (e.g. 

; IThomas et~aH |2005|; iBernardi et all 
2006; iThomas et aJj|201Ch . Most intere- 



Trager et all [2000b 



2006; Clemens et al 



stingly, the [a/Fe] ratio participates in the E-E dichotomy, 
i.e. elliptical galaxies with low [a/Fe] ratios have core-less 
central profiles, while [q/Fe]-enh anced galaxies with sh ort 
formation time-scales have cores (|Kormendv et al.ll2009h . 

Individual element abundance ratios, in addition to 
a/Fe, can further disentangle the formation of different 
stellar populations. Since the individual elements are pro- 
duced in different stellar e volutionary phases they trace 
varying form a tion histories, fedvards son et al.l (|l993T ) and 
iBensbv eT al. (2010) derive as many as 12 different element 
abundance ratios to distinguish between the different forma- 
tion histories of the stellar populations in the Milky Way. 

A number of studies have derived individu- 
al element abundance ratios for early-type galaxies 



dSanchez-Blazauez et all 120031, l2006al; 


Clemens et al.l:2006; 


Kelson et al.l 20061; Graves et al.l 2007 


; Graves & Schiavonl 


2008; Smith et alj|2009; Price et al. 201lf), in several cases 



for fairly small samples. Different methods are applied 
in these studies, but they are all based on absorption 
line indices. The results are dependent on the method 
applied and the sample used. The aim of this work is to 
simultaneously derive all element abundance ratios allowed 
by the sensitivity of the adopted absorption line indices. 
The maximum amount of information is extracted from the 
indices to reliably derive the abundance ratios and state of 
the art models of stellar populations of absorption indices 
are utilised to obtain as accurate results as possible. 

To fully interpret observed element abundance ra- 
tio trends, stellar nucleosynthesis needs to be understood. 
iPipino et all (|2010T ) find up-to-date models of chemical evo- 
lution to struggle in simultaneously reproducing observed 
abundance ratios for Carbon and Nitrogen from the unresol- 
ved stellar populations of early-type galaxies. The a-element 
Ca is a puzzle as it is has been found to tr ace Fe instead 



of other a-elements for early-type galaxies (jCen arro et al 



Saglia et al 



20031 ; IThomas etakl l2003bl ; ISmith et all [200' 
20021 '). This has been inte rpreted as Ca being contributed 
by SN Ia as well as SNII (jThomas. Johansson fc MarastorJ 
l2011al ). Thus element abundance ratios of unresolved stellar 
populations are also useful for constraining stellar nucleo- 
synthesis. 

We present a technique for deriving a wide range of 
element abundance ratios for unresolved stellar populati- 
ons, including [O/Fe] (representing [a/Fe]), [C/Fe], [N/Fe], 
[Mg/Fe], [Ca/Fe] and [Ti/Fe]. The method is based on new 
flux-calibrated stella r population models of absorption line 
indices presented in IThomas et ail (|2011bl ). We analyse a 
sample of 3802 SDSS early-type galaxies for which we inve- 



stigate element ratio scaling relations with velocity disper- 
sion. 

The paper is organised as follows. The SDSS early-type 
galaxy sample used is presented in Section [2] and the techni- 
que for deriving the element abundance ratios is described in 
Section [4] The results of derived element abundance ratios 
for the data sample are presented in Section \S\ and further 
discussed and compared with the literature in Section [6] 
Conclusions are given in Section [7] 



2 THE DATA SAMPLE 

The selected sample is part of the MOSES catalogue (Mor- 
phologically Sel ected Early-types in SD SS). This is descri- 
bed in detail in ISchawinski et al.l (|2007l ) and IThomas et ail 
(|2010T ) (T10) and only a brief description is given here. The 
MOSES sample consists of 48 023 galaxies from the SDSS 
Data Release 4, selected to have a magnitude r < 16.8 in 
the redshift range 0.05 < z < 0.10. The selection criteria 
ensured a reliable visual inspection of galaxy morphology 
and the full sample was divided into late-types (31 521) and 
early- types (16 502) through a purely visual classification 
scheme. 

As comprehensively described in T10 the 25 standard 
Lick absorption line indices were measured on the galaxy 
spectra downgraded from the SDSS spectral resolution to 
the Lick/IDS resolution and corrected for emission line fil l 
in using G ANDA LF/PPXF (|Cappellari fc Emselleml |2004| ; 
ISarzi et al.ll2006T l. The stellar kinematics and best-fit stellar 
templates from PPXF were used in T10 for correcting the 
measured indices for stellar velocity dispersion broadening 
effects. The velocity dispersion measurements derived and 
published in T10 are used in Section [5] as proxies for galaxy 
mass. 

The visual classification does not bias against star for- 
ming galaxies and the sub-sample of early-type galaxies will 
therefore include galaxies with blue colours, having possible 
on-going star formation or recent star formation. These ga- 
laxies can have possible emission line contamination in the 
absorption features, but the use of GANDALF gives reliable 
Lick index measurements for the full sample. 

The final sample (3802 objects) used in this work was 
selected according to T10. This is a sub-sample of early-type 
classified galaxies in the narrow redshift range 0.05 < z < 
0.06, chosen to minimise evolutionary effects present when 
using the full range of redshifts, affecting the derivation of 
stellar population parameters. 



3 THE TMJ MODELS 



In IThomas. Maraston fc Johanssonl (|2011br i (TMJ) we pre- 
sent new stellar population models of Lick absorption-line 
indices with variable element abundance ratios. The model 
is an extension of the lThomas et aL | (|2003al . l2004h (TMB/K) 
model, which is base d on the e v olutionary s tellar populati- 
on synthesis code of iMarastonl (|l998L 120051) . A calibration 
on gala ctic globular clusters was performed in lThomas et al.l 
(|2011aT ) for the range of element abundance ratios covered. 
For basic information on the model we refer the reader to 
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iThomas etafl |2003al . l2004h and TMJ. Here we provide a 
brief summary of the main features of our new models. 



3.1 New features 

The key novelty compared to the TMB model is that the 
TMJ model is flux-calibrated, hence not tied anymore to the 
Lick/IDS system. This is because the new models are based 
on our calib rations of absorption- line indices with stellar 
parameters dJohansson et alll2010l ) derived from the flux- 
calibra ted stellar library MILES l|Sanchez-Blazguez et alj 
2006b). The MILES library consists of spectra covering the 
wavelength range 3500-7400 A and with a spectr al reso- 
lution of 2.54 A (as revised by iBeifiori et alj|201ll ) of 985 
stars selected to produce a sample with extensive stellar pa- 
rameter coverage. Most importantly it has been carefully 
flux-calibrated, making standard star-derived offsets unne- 
cessary. The data release now provides two model versi- 
ons partially based on different stellar evolutionary tracks, 
ICassisi et all (|l997T ) and Padova (|Girardi et alj2000h at high 
metallicities. In the present study we use the former version, 
since these tracks were used in the TMB/K model adopted 
in TIP, it wa s calibr ated on galactic globular clusters in 
I Thomas et alj (|2011al ) and our Padova model version only 
cover ages down to 2.8 Gyr. 

Finally, we calculate models in which we selectively en- 
hance, by 0.3 dex, each of the elements C, N, Na, Mg, Si, 
Ca, Ti, and Cr in turn, released with TMJ. Further exten- 
ded model calculations with arbitrary enhancement factors, 
hence a finer grid in element ratios, are used in the present 
work. 



3.2 Element abundance ratios 



As in lThomas et al.1 l|2003af ) and IThomas et all ( |2011bl ) we 
keep the total metallicity fixed while varying the element 
abundance ratio [a/Fe]. Briefly, Fe-peak elements (Fe and 
Cr) are locked together (depressed group) and the rest of 
the considered elements (C, N, O, Na, Mg, Si, Ca and Ti) 
are locked to the a-elements (enhanced group). Solar va- 
lues of element ratios a re known from measurements of the 
individual abundances (|Grevesse. Noels fc Sauva3ll996r ). A 
non-solar [a/Fe] ratio is computed by simultaneously chan- 
ging the abundances of elements in the enhanced and de- 
pressed group to counter-balance a change of total metal- 
licity. A non-solar [a/Fe] ratio mainly means a change in 
the abundances for the Fe-peak elements, since a-elements 
and in particular O xygen dominate tota l me tallicity in the 
Sun, as discussed in lTrager et al.l l|2000al ) and IThomas et all 
(|2003ah . 

We keep the total metallicity fixed also when we pro- 
duce models with enhancements of the individual elements 
C, N, Mg, Ca and Ti. Starting from an [a/Fe] ratio the 
models are perturbed by enhancing/depressing the indivi- 
dual element E. Hence we actually consider elements rati- 
os of given element over a, e.g. a N-enhanced model ha- 
ve varying [N/a] ratios. The [E/Fe] ratio is then calculated 
with [E/Fe] = [E/a] + [a/Fe], where [a/Fe] is the starting ele- 
ment ratio. The element E is detached from the rest of the 
elements in the enhanced group, while the ratio between 
elements in the enhanced and depressed group is locked. 



Thus when enhancing/depressing the individual element E 
the total metallicity is conserved and the abundance ratios 
between the rest of the elements remain unchanged. Vary- 
ing [E/Fe] mainly means a change in the abundance of the 
element E, since the rest of the elements, locked together, 
dominate total metallicity. 

3.3 Re-calculation of [Fe/H] 

In addition to the pa rameters discussed above we de - 
rive [Fe/H] Follow ing iTantalo. Chiosi fc Bressanl (|l998h . 
iTrager et all (|2000ah and lThomas et all (|2003ah we have the 
relationship between iron abundance [Fe/H] and total metal- 
licity [Z/H] when only considering the element ratio [a/Fe] 

[Fe/H] = [Z/H] + A[a/Fe] = [Z/H] + A[0/Fe] (1) 

if we assume that a/Fe reflect O/Fe (see Section l"3.4[) . Hence 
Eq. [T] remains valid if the analysis is restricted to [a/Fe] 
([O/Fe]). Eq. [T]instead needs to be revised when considering 
more element ratios besides [a/Fe] ([O/Fe]). Modifying the 
individual element ratios with [O/Fe] as a starting point as 
in this work (see Section [3.2[) . perturbations to Eq. [T] arise 
in the form [E/O] for element E. We rewrite the general 
relationship of Eq. fTJto 

[Fe/H] = [Z/H] + A[0/Fe] + ^ B x [E x fO] (2) 

X 

for x number of individual elements E^ with corresponding 
coefficients B^. Changing an individual element abundance 
while keeping total metallicity and the ratio between the 
rest of the elements fixed, Eq. [5] becomes 

A\Fe/H] = B x A[E x /0] = B X (A[E X /H] - A[0/H]) (3) 
leading to 

B x = (4) 

A[E x /H]-A[0/H] {) 

Following IThomas et al.1 l|2003al '). when varying the [a/Fe] 
ratio only, total metallicity is kept fixed while enhancing 
the a-elements through 



f a X+ + f Fe x- = x+ + x~ 



(5) 



where X + and X - are the mass fractions of the enhanced 
and depressed groups (see Section I3.2|l , respectively, chan- 
ged by the factors f a and fFe- Hence for an increase in 
the [a/Fe]-ratio a higher abundance of the a-elements is 
counter-balanced by a decrease in the abundance of Fe-like 
elements to keep total metallicity fixed. Further, changing 
an individual element abundance keeping total metallicity 
and the ratio between the rest of the elements fixed, Eq. [5] 
is extended to 

f E (f a X E ) + MfFeX- + fa(X + — X E )) = 

f a X+ + f Fe X- — X + + X~ (6) 

where the individual element abundance is now multiplied 
by the factor fg and the rest of the elements by fo . Hence if 
the abundance of element E is changed by the factor fe, the 
abundances of the rest of the elements are changed by the 
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Figure 1. The response of the 25 Lick indices to changes in individual element abundances for a 12 Gyr solar metallicity modelled 
stellar population. The response is defined as the fractional index change s for enhancements of a specific element abundance by a factor 
two and normalised to the typical observational index errors presented in [Johansso n et all d2010T) . For each element we display the index 
responses on an individual scale (left-hand panel) as well as those referred to a fixed scale (right-hand panel, see text for details). 



factor fo to counter-balance a change in total metallicity. 
Also, since all other elements besides E are changed by the 
same factor the ratios between these elements remain con- 
stant. This implies secondary abundance changes due to f_B 
and fo besides that of i a and ffe- The logarithmically solar 
scaled iron abundance is defined as 



[Fe/H] old = log 



X 1 



X H 



log 



(7) 



where X Fe and X H are the mass fractions of Fe and H, 
while Xq e and Xq are the corresponding solar values. If the 
mass fraction is changed by the factor fo the iron abundance 
becomes 



[Fe/ff]„ ew =log(^)+logC^ 



log 



Xr. 



A 



(8) 



Eq. tag give A[Fe/H]=[Fe/HW[Fe/H] o;d =log(/ ). In the 
same way we get A[0/H]=log(/o) and A[Ex/H]=log(/E) 
such that Eq. [4] becomes 



B x 



log(/c 



log(/ B ) - log(/o) 



(9) 



Similarly to Eq.[7][S]we can enhance the element ratio [E/0] 
from 



\E/O\oid = log 



to 



[E/0]new = log 



+ log 

+ log 
log 



X° 



hi 

fa 

x° ) 



+ 



log 



*8 



(ii) 



Table 1. Solar abundance fractions (Xg) derived from 
iGrevesse et al.l |l996), abundance enhancement factors (f^; and 
fo) and final element ratio coefficients (B x ) for the various ele- 
ments considered. 





x 


fe 


fo 


B,, 


c 


0.172 


1.693 


0.847 


-0.24 


N 


0.053 


1.894 


0.947 


-0.079 


Mg 


0.038 


1.923 


0.962 


-0.056 


Ca 


0.003 


1.994 


0.997 


-0.0043 


Ti 


0.0002 


1.9996 


0.9998 


-0.00029 



For a new abundance ratio [E/O] ne iu=0.3 starting from solar 
values [E/O] oid =0.0 Eq.QDand Eq.lTTIgive 



log ^ 



0.3 



(12) 



Following iThomas et "all (|2003al ) with the difference of ha- 
ving C locked to the enhanced group we re-derive fFe 
and f a . Adopting the solar abundance frac t ions X + =0.91 
and X~ =0.079 derived from IGrevesse et all l|l996t) we find 
fFe=0.521 and f a = 1.042, resultin g in A=0.94. Individua l ab- 
undance fractions derived from IGrevesse et al.l (Il996t ) are 
presented in Table [T] along with the corresponding values 
for fs and fo derived using Eq.[6]and Eq. 1121 Eq.[9]then gi- 
ves the final coefficients B^, also presented in Table [1] The 
final relationship becomes 

[Fe/H] = [Z/H]-0.94[O/Fe\-0.2A[C/O]-0.079[N/O]- 
0.056[Mg/O] - 0.0043[C*a/O] - 0.00029[Ti/O] (13) 
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3.4 Index Responses 

Fig.[T]shows the response of the 25 Lick indices to individual 
element abundance changes for a 12 Gyr, solar metallicity 
stellar population. The fractional index change is calculated 
for an enhancement of the respective element by a factor two 
normalised to the ty pical observationa l meas urement error 
for MILES stars from [Johansson et al, (2010). The scale on 
the x-axis in the left hand panel of Fig.[T]is different for each 
element. Hence the left-hand panel should be read vertically 
demonstrating the most sensitive indices for the individual 
elements. The scale on the x-axis is instead kept fixed for 
all elements in the right hand panel. This panel should then 
instead be read horizontally to identify easily those elements 
that are best traced by a specific index for the current set 
of models. It can be seen that the elements C, N, Na, Mg, 
Ca, Ti, and Fe are best accessible. 

The abundance of nitrogen is obtained from the CN in- 
dices that are also highly sensitive to C abundance. However, 
this degeneracy can be easily broken through other C sensi- 
tive indices such as C24668 and Mgi. The Mg indices Mgi, 
Mg2, and Mgb are very sensitive to Mg abundance. Note, 
however, th at all three additionally anti-correlate wi th Fc 
abundance (|Trager et al ] |2000al ; iThomas et al.|[2003al ). Ca 
can be measured well from Ca4227, except that this parti- 
cular index is quite weak and requires good data quality. Na 
abundance can be derived quite easily from NaD in princi- 
ple. However, in practise this is problematic as the stellar 
component of this absorption feature is highly con tamina- 
ted by interstellar absorption (|Thomas et al.ll2003al ). Iron is 
well sampled through the Fe indices. 

There are two among the Fe indices, however, that are 
also sensitive to Ti abundance besides Fe. These are Fe4531 
and Fe5015. They offer the opportunity to estimate also Ti 
abundance. We will only use Fe4531, as Fe5015 is contami- 
nated by a non-negligible Mg sensitivity besides Fe, which 
weakens its usefulness for Ti abundance determinations. 

The remaining three elements O, Si and Cr cannot ea- 
sily be measu red through the availab le indices. As discussed 
extensively in I Thomas et all (1200381 ). however, oxygen has a 
special role. O is by far the most abundant metal and clearly 
dominates the mass budget of 'total metallicity'. Moreover, 
the a/Fe ratio is actually characterised by a depression in 
Fe abundance relative to all light elements (not only the a 
elements), hence a/Fe reflects the ratio between total me- 
tallicity to iron ratio rather than a element abundance to 
iron. As total metallicity is driven by oxygen abundance, the 
a/Fe can be most adequately interpreted as O/Fe ratio. We 
therefore re-name the parameter a/Fe to O/Fe under the as- 
sumption that this ratio provides an indirect measurement 
of oxygen abundance. 



4 DERIVATION OF STELLAR POPULATION 
PARAMETERS 

The aim of this work is to develop a method for deriving 
a comprehensive set of element ratios for unresolved stellar 
populations to assess the early-type galaxy sample descri- 
bed in Section [5] This continues the work of T10 who de- 
termined the stellar population parameters Age, [Z/H] and 
[a/Fe]. Including these parameters we extend the work by 



simultaneously deriving the individual element abundance 
ratios [C/Fe], [N/Fe], [Mg/Fe], [Ca/Fe] and [Ti/Fe] using 
the stellar population models described in Section We in- 
clude the element Ti, one of the heaviest a-elements before 
the Fe-group. In the following sections we introduce the me- 
thod. 

4.1 Method 

First we determine the traditional light-averaged stellar po- 
pulation parameters age, total metallicity, and [O/Fe] (re- 
presenting [a/Fe], see Section 13. 2p from a base set of in- 
dice s, using a y 2 minimiz ation routine (described in T10 
andrTh onias et all l|2011ah ). This base set of indices inclu- 
des Mg b, the Balmer indices Hsa, H^f and H/3 and the iron 
indices Fe4383, Fe5270, Fe5335, Fe5406, chosen as they are 
well calibrated with galactic globula r clusters without ind ivi- 
dual element abundance variations (|Thomas et al.ll2011al lbT) . 
The models used at this point have step sizes of 0.024 
dex for log(age) and 0.05 dex for both [Z/H] and [O/Fe], 
over the ranges 0<age<20 Gyr, -2.25<[Z/H]<0.7 and - 
0.3<[O/Fe]<0.5. Only indices that are sensitive to these 
three parameters are included in the base set. Mgb is inclu- 
ded in the base set of indices since it is usefu l in combinati- 
on wit h Iron indices to constrain [O/Fe] (e.g. iThomas et all 
l2003al ). 

In the subsequent steps we add in turn particular sets 
of indices that are sensitive to the element the abundance 
of which we want to determine (see step 2-6, Fig. [2}. In 
each step we re-run the x 2 fitting code with a new set of 
models to derive the abundance of this element. This new 
set of models is a perturbation to the solution found for the 
base set of indices. It is constructed by keeping the stellar 
population parameters age, metallicity, and O/Fe fixed and 
by modifying the element abundance of the element under 
consideration by ±1 dex in steps of 0.05 dex around the base 
value. 

A new best fit model is obtained from the resulting x 2 
distribution. Then we move on to the next element. Due to 
the fact that several indices respond to the same elements 
(see Fig. [1} an iterative method is needed to simultaneous- 
ly derive the individual element abundance ratios, illustra- 
ted in Fig. [2] The derivation of individual abundance ratios 
(Step 2-6) is iterated until the abundance ratios remain un- 
changed within the model grid step of 0.05 dex. A fast con- 
vergence of 3-4 iterations is generally found for this inner 
loop and a maximum of 5 iterations is set as an upper limit. 

At the end of the sequence we re-determine the overall 
X 2 and re-derive the base parameters age, metallicity, and 
O/Fe for the new set of element ratios. At this final step we 
use models with the base parameters age, total metallicity, 
and O /Fe ratio perturbed around the previously derived va- 
lues by Alog{age)/A[Z/H}/A[0/Fe] = ±0.1 dex and with 
step sizes of 0.02 dex for log(age) and [Z/H] and 0.05 dex for 
[O/Fe]. A bigger step size for [O/Fe] improves upon speed 
for the routine and it was found that it did not affect the fi- 
nal results. All indices considered in the inner loop (step 2-6) 
are used at the final step together with the base set of in- 
dices. Then we go back to the second step and use the new 
base parameters to derive individual element abundances. 
This outer loop is iterated until the final x 2 stops improving 
by more than 1 per cent. The method converges relatively 
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Figure 2. Method to derive element abundance ratios. A x 2 -minimisation routine is used at all steps to find the best fit model. Step 1 
computes an initial fit for the base parameters age, [Z/H] and [O/Fe] using a base set of indices sensitive to these parameters only. The 
base parameters are fixed when deriving the individual element abundance ratios (inner loop, step 2-6) by adding indices sensitive to 
the element considered in each step to the base set of indices. At convergence, i.e. the element abundance ratios remain unchanged, the 
routine exit the inner loop and the base parameters are re-derived for the computed set of abundance ratios using all indices. This outer 
loop is iterated until the \ 2 a t step 7 stops improving by less than 1%. Models with varying abundance ratios of the element considered 
are produced in step 2-6, while models perturbed around the previously derived base parameters are produced in step 7. The indices 
used at each step are presented in each box and the typical number of iterations are indicated as convergence along with the maximum 
number iterations allowed for both the inner and outer loop. 



fast, again generally requiring 3-4 iterations and 5 iterations 
is set as an upper limit. 

In more detail, the sequence of elements is as follows. 
The first element in the loop is carbon, for which we use 
the indices CN1, CN2, Ca4227, G4300, H 7 a, H 7 f, C24668, 
Mgi, and Mg2 on top of the base set. Next we drop the- 
se C-sensitive indices and proceed deriving N abundance, 
for which we use the N-sensitive indices CNi, CN2, and 
Ca4227. Then we move on to Mgi and Mg2 for Mg ab- 
undance, Ca4227 for Ca, and finally Fe4531 for the element 
Ti. 

At step 1 and 7 we allow the procedure to discard indi- 
ces with a bad x 2 > explained in detailed in T10. Briefly, the 
probability distribution around the ^-minimum is compu- 
ted with the incomplete F-function for the model grid used 
at each step. This distribution gives the probability P that 
the true model has a x 2 lower than the value obtained in the 
fit. If P ^ 0.999 the fit is considered to be unacceptable and 
the index with the largest \ 2 i s discarded. This procedure 
is repeated until P > 0.999. In 5.2% of the cases one index 
was discarded, while more than one index was discarded for 
1.6% of the 3802 objects. Compared to T10 the percentage 
where at least one index is discarded has been reduced from 
>30% to 5.8%. This is partly due to the fact that some indi- 



ces used in T10 have not been considered here, e.g. Ca4455 
and NaD that were the most frequently discarded indices in 
T10. The inclusion of individual element abundance ratios 
improve the fit for several indices and is thus also responsible 
for the lower number of discarded indices. Table [2] gives the 
discard percentage for all Lick indices. G4300 is the index 
most frequently removed, as it was discarded in 3.8% of the 
cases, while Fe5270 reached 1.1% and the rest of the indices 
had a discard percentage less than 1.0%. 

4.2 Errors 

Errors on the parameters are estimated by taking the 
FWHM of the probability distributions (see previous sec- 
tion) and converting to lcr errors using <r=FWHM/2.355. 
Alog(age), A[Z/H] and A [O/Fe] are derived at step 
1 (see previous section) since the models grids are too 
narrow at step 7 to reliably estimate the errors. Sin- 
ce the individual abundance ratios are derived as per- 
turbations to the [O/Fe] ratio, errors on these parame- 
ters are estimated by in quadrature adding the errors 
on the derived perturbations A[E/Fe] p to A[0/Fe], i.e. 
A[E/Fe\ 2 =A[0/Fe] 2 +A[E/Fe] 2 p . Fig.[3]shows the norma- 
lised probability distributions for 100 randomly selected ob- 
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Figure 3. Column 1 and 3: Probability distributions for 100 randomly selected objects (black gaussians) and 4 highlighted objects 
with varying parameter errors (coloured gaussians). Column 2 and 4: Distributions of parameter errors, which are estimated from the 
probability distributions (see text for more details), for the the full sample. Black dashed line shows the average errors, while colour 
dashed lines shows the errors for the highlighted probability distributions in column 1 and 3 with corresponding colours. 



jects (Eiffel tower plots in 1st and 3rd column, black gaus- 
sians) of our sample (see Section [2) together with the dis- 
tributions of the estimated errors (2nd and 4th column). 
Four objects are highlighted demonstrating the relationship 
between the sizes of the errors and the widths of the proba- 
bility functions. The dashed vertical lines in column 2 and 
4 show the error estimates of these objects and are coloured 
according to the highlighted probability distributions in co- 
lumn 1 and 2. Black dashed vertical lines in column 2 and 4 
are average errors. The probability distributions are moved 
to a common position by shifting them with the parameter 
value at peak probability (P roM ). In general we find very 
well defined one-peaked gaussians. For increasing errors the 
probability distributions clearly increase in width. 

Fig' |3] shows the relationship between the estimated er- 
rors for the base parameters (age, [Z/H] and [O/Fe]) and the 
signal-to-noise (S /N) ratio of the galaxy spectra for the 3802 
early-type galaxies studied in this work (contours). The left 
hand and right hand panels show the relationship for S/N 
in the g-band and r-band, respectively. The S/N ratios are 



the median values of S/N per pixel in each band as given 
by the SDSS spectroscopic fits-headers. Orange solid lines 
are mean errors in bins of S /N with a width of 5 dex. Below 
a S/N of ~30 and ~45 for the g-band and r-band, respec- 
tively, the error-S/N relationships are close to linear. Thus 
the green dashed lines show linear relationships in these re- 
gimes on the form error=axS/N+b with the values of the 
fit parameters a and b given in Table [3] 

4.3 Globular cluster calibration 

The method describe d in the previous section is used in 
iThomas et all (|2011al ) to obtain estimates of stellar popu- 
lation parameters and element ratios from integrated light 
spectroscopy of galactic globular clusters. We show that the 
model fits to a number of indices improve considerably when 
various variable element ratios are considered. Our derived 
ages, metallicities and abundance ratios agree generally very 
well with the literature values from photometry and high- 
resolution spectroscopy of individual stars. 
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Figure 4. Relationship between the estimated errors of the base parameters age, [Z/H] and [O/Fe] with S/N ratio in the g-band (left 
hand panel) and r-band (right hand panel) for the 3802 early- type galaxies studied (contours). Orange solid lines are mean errors in bins 
of S/N with a width of 5 dex and green dashed lines show linear relationships below a S/N of ~30 and ~45 for the g-band and r-band, 
respectively, on the form error=axS/N+b. The fit parameters a and b are given in Table|3] 



Table 2. Discard percentage for all Lick indices 



Index 


Discard % 


HiA 


< 1.0 


H 5jP 


< 1.0 


CNi 


< 1.0 


CN 2 


< 1.0 


Ca4277 


< 1.0 


G4300 


3.8 


H 7 a 


< 1.0 


H 7 F 


< 1.0 


Fe4383 


< 1.0 


Ca4455 


not used 


Fe4531 


< 1.0 


C 2 4668 


< 1.0 


H/3 


< 1.0 


Fe5015 


not used 


Mgi 


< 1.0 


Mg 2 


< 1.0 


Mgb 


< 1.0 


Fe5270 


1.1 


Fe5335 


< 1.0 


Fe5406 


< 1.0 


Fe5709 


not used 


Fe5782 


not used 


NaD 


not used 


TiOi 


not used 


TiO a 


not used 



4.4 Comparison with methods in the literature 

iTrager et all (|2000al ) introduced models with varying ele- 
ment abundance ratios by using t he single stellar population 
models (SSPs) of lWorthevI (I1994T) together with the Lick in- 
dex response functions of lTripicco fc Belli (|1995T ). They inve- 
stigate the nature of the elements O and C by using different 
model treatments of these elements, i.e. locked to the enhan- 
ced group (see Section 13.2ft . locked the depressed group or 
kept fixed at solar valu es. panchez-Blazquez et al.l (|2006al ) 
follow t he technique of ITrager et alj (|2000al ) , but use the 
SSPs of IVazdekis et all i|2010r i" and modelled indices with 



Table 3. Linear error-S/N relationships below S/N~30 and 
S/N~45 for the g-band and r-band, respectively, on the form 
error=axS/N+b. 



Band 


Parameter 


a 


b 


g 


A log (age) 


-0.0077 


0.38 


g 


A[Z/H] 


-0.0075 


0.35 


g 


A [O/Fe] 


-0.0043 


0.18 


r 


A log (age) 


-0.0062 


0.43 


r 


A[Z/H] 


-0.0047 


0.34 


r 


A [O/Fe] 


-0.0038 


0.22 



different treatments for C, N and Mg to find the description 
that in general best resembles the overall behaviour of the 
dat a considered. 

ISanchez-Blazquez et alj (|200Sf ) indirectly predict diffe- 
rences in the element abundances of C and N. They use 
absorption line indices known to be sensitive to variations 
of the abundances of these elements. If these indices show 
differences, while indices insensitive to abundance variations 
of the elements considered show no difference the stellar po- 
pulations are believed to have differ ent element abundanc e 
ratios. They also consider the SSPs of lVazdekis et all |l996), 
but without elemen t abun dance ratios. 

IClemens et all ll2006fl treat [C/H] as a separate va- 
riable and conseque ntly fit modelled indices (SSPs from 
lAnnibali et al.l 2007) wit h the free paramete rs age, [Z/H], 
[a/Fe] and [C/H] to data. iKelson et all (|2006h consider mo- 
delled indices (TMB/K) with the free parameters age, [Z/H], 
[a/Fe] and [a/C], [«/N] and [a/Ca], but not all parameters 
simultaneously. 

The first authors to develop a method that simulta- 
neously consider a full grid of element abundance r atios were 
iGraves fc Schiavonl (l2008fl . based on the models of lSchiavonl 
( 20071 ). They start by deriving a fiducial Age and [Fe/H] 
using the index H/3 in combination with <Fe>. H/3 is then 
exchanged with an index sensitive to a specific element. This 
element is then either enhanced or depressed until the fidu- 
cial Age and [Fe/H] is matched. 
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Figure 5. Comparison between the derived parameters of this work (y-axes) and T10 (x-axes) for age (left hand panel), total metallicity 
(middle panel) and [O/Fc] (right hand panel). For the latter we compare to the [a/Fe] ratios from T10 as it corresponds to [O/Fe] of this 
work (see Section |3,4|I , Contours show the density of data points and solid lines are least-square fits. These results are derived using the 
method presented in this work (Jll setup, see Section[4]l. For the TMJ models (see Section[3} we have also derived the stellar population 
parameters with the same setup as used in T10. Least-square fit representing the relationship between these results and T10 are shown 
as dot-dashed lines (T10 setup). Fit parameters are given by the labels and dashed lines are 1-to-l relationships 



The method developed in this work is based on a diffe- 
rent philosophy. In our method a \ 2 routine finds the model 
with an enhancement of a specific element that best fits 
the data keeping a set of base (fiducial) parameters (Age, 
[Z/H] and [O/Fe]) fixed. Initially deriving [O/Fe] together 
with age and [Z/H] gives more freedom to the choice of age- 
sensitive indices. Thus we can use the Fe-sensitive higher 
order Balmer indices (H<5a, H<5f, H7A and ~Hjf) besides H/3 
for initially constraining age. A major difference is that we 
allow for iteration between the base model and the soluti- 
on with varying element ratios. For the sample described in 
Section [2] the average difference between the initially deri- 
ved base parameters and the final values are 0.11 dex for 
log(age), 0.09 dex for [Z/H] and 0.04 dex for [O/Fe]. 

For each element we use all indices that are sen- 
sitive to its abundance (see Fig. [T] Section 13. 4p . while 
iGraves k, Schiavonl l|2008l ) use the most sensitive index only 
for each element. Using all indices extracts all information 
available and protects against anomalies in individual indi- 
ces, i.e. noise and emission line fill affecting the absorpti- 
on features. Still, more sensitive indices will have a greater 
weight than less sensitive indices. [Graves fc Schiavonl (2008) 
set a fixed value of [O/Fe] and simultaneously derive [C/Fe], 
[N/Fe], [Mg/Fe] and [Ca/Fe]. We extend this and further in- 
clude [Ti/Fe], while we also trace [O/Fe] inferred from [a/Fe] 
(see Section [3. 4p . 



5 RESULTS 

In the following we present stellar population parameters 
and element abundance ratios as functions of the stellar velo- 
city dispersion measurements derived in T10 (see Section[2]). 
Since this is a continuation of the work of T10 we compa- 
re the results of the base parameters (age, Z/H and O/Fe) 
with their results, while the new element ratios considered 
are compared with the full literature. 

5.1 A direct comparison with T10 

In Fig.[S]we compare the ages, [O/Fe] ratios and total metal- 
licities derived using the method presented in Section [4] (Jll 
setup) to the corresponding parameters from T10 (where 



[a/Fe] corresponds to [O/Fe] of this work, see Section [3.4p . 
Contours show the density of data points and dashed li- 
nes are 1-to-l relationships. Solid lines are least-square fits 
with fit parameters given by the labels. We find a very good 
agreement for the derived ages and generally good agree- 
ments for total metallicity and [O/Fe]. For the latter two 
this work shows overall higher values and the differences 
increase towards lower parameter values. Still, the average 
differences are small within the parameter range containing 
the majority of data points (<0.2 dex for [Z/H] and <0.1 
dex for [O/Fe]) . 

Using the TMJ models (see Section[3]) we have also deri- 
ved the stellar population parameters age, [a/Fe] and total 
metallicity using the same setup as used in T10 (TMJ25 
setup), i.e. using all 25 Lick indices and not considering ele- 
ment ratios beyond [a/Fe]. With these results we can disen- 
tangle the effect on the derived parameters from using both 
new models and accounting for individual element ratios. A 
direct comparison between the TMJ25 and T10 is also in- 
cluded in Fig. [5] (least-square fit dash-dotted lines). Again 
we find a good agreement between the derived ages. In the 
old age regime the TMJ25 ages are somewhat higher (~0.1 
dex). Hence in this regime the TMJ models seem to result in 
slightly older ages compared to the TMB/K models, while 
accounting for element ratios results in younger ages. The 
TMJ25 and T10 metallicities are very similar within the pa- 
rameter range containing the majority of data points (con- 
tours). Hence, the metallicities do not change considerably 
when using the new models. Instead, the metallicities are 
slightly lower if individual element ratios are not taken into 
account. 

For [a/Fe] TMJ25 produces overall higher element ra- 
tios by ~0.1 dex compared to T10. In this case the choice of 
models matters as well as accounting for individual element 
ratios. In the right hand panel we have also included a com- 
parison between the Jll [Mg/Fe] ratios and the [a/Fe] ratios 
from T10 (least-square fit blue line). In fact the [Mg/Fe] ra- 
tios agree very well with the TMJ25 [a/Fe] ratios as the blue 
and dash-dotted are very similar. Hence the [a/Fe] ratios de- 
rived without taking individual element ratios into account 
reflect [Mg/Fe] rather than [O/Fe]. This is probably due to 
the fact that the Lick indices are more sensitive to variations 
in the abundance of Mg than O (see Section \3A\ . 
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Figure 6. Left hand panel: shows the relationship between the derived light-average ages and velocity dispersion. The dashed orange 
line indicate the separation of an old r ed sequence populat ion (orange contours) from a rejuvenated population (cyan contours) with 
light-average ages smaller than 2.5 Gyr llThomas et aLlfeoiOf) . The fraction of rejuvenated galaxies is indicated below the separation line. 
The whole sample is shown as grey-scaled filled contours. The orange solid line is a least-square fit to the red sequence population and 
the parameters of the fit are given at the top of the panel. The dashed black line is the analogous fit from T10 for comparison. Median 
lcr-errors are shown in the lower right corner. Right hand panel: shows the distribution of deviation in age from the least-square fit to 
the red sequence population with the same colour coding as in the left hand panel. The standard deviation of the gaussian fitted to the 
distribution of the red sequence population (orange line) is indicated in the upper right corner. 



The TMJ25 setup consider all Lick indices and the Jll 
setup a selection of 18 indices best calibrated with galactic 
globular clusters (see Section[4]). The difference seen between 
these two setups may therefore arise from the choice of in- 
dices instead of accounting for element ratios. To evaluate 
this we also derived ages, [a/Fe] ratios and total metallici- 
ties using the TMJ models without element ratios, but for 
the 18 indices considered in this work. These result are very 
similar to that of TMJ25, such that the different choices of 
indices do not matter significantly. 



5.2 Ages 

The relationship between luminosity-weighted age and velo- 
city dispersion is presented in the left hand panel of Fig [6] 
The full sample is shown with grey-scaled, filled contours. 
We reproduce the result from T10 of having a bimodal dis- 
tribution of ages in analogy to a red sequence populati- 
on (log(age)>0.4, orange contours) and a rejuvenated blue 
cloud population (log(age)<0.4, cyan contours) of low mass 
early-type galaxies, identified thanks to the purely visual 
classification of the MOSES sample (see Section [2j. The 
fraction of rejuvenated galaxies is 13.70%, close the corre- 
sponding fraction of 10.15% found in T10. The rejuvenated 
population is further discussed in the following sections and 
in more detail in Section 16.11 Fig. [6] also shows a cr-clipped 
(3cr limit) linear least-square fit to the red sequence popu- 
lation (orange solid line), along with the analogous from 
T10 for comparison (dashed black line). We find very simi- 
lar trends, but with a slightly steeper slope, as in T10, i.e. 
increasing ages with increasing velocity dispersion. 

The distribution around the fit to the red sequence po- 
pulation is shown in the right hand panel of Fig. [6] for the 
whole sample (grey histogram), red sequence population and 
rejuvenated population (black open histograms). Gaussians 
fitted to the rejuvenated and red sequence distributions are 



shown with the same colour coding as in the left hand pa- 
nel, with the standard deviation of the latter given by the 
label. The rejuvenated population shows a distinct peak, 
offset from the red sequence population to younger ages by 
~0.55 dex. 



5.3 [Z/H] 

The left hand panel of Fig. [7| shows the relationship bet- 
ween total metallicity and velocity dispersion, with the sa- 
me colour coding as in Fig. [S] The least-square fit to the red 
sequence population (orange line) indicate a strong correla- 
tion between total metallicity and velocity dispersion. This 
is in agreement with T10, but the slope of the fit derived 
here is smaller than the corresponding fit derived in T10 
(black dashed line). In Section \5. II we find that not accoun- 
ting for element ratios produce the lower metallicities found 
in T10 compared to this work. This difference increases to- 
wards lower metallicities and produce the flatter slope of the 
[Z/H]-<7 relation. 

The distribution around the fit to the red sequence po- 
pulation is shown in the right hand panel of Fig. [7) The reju- 
venated population (cyan gaussian) shows an offset towards 
higher metallicities by ~0.1 dex compared to the red se- 
quence population (orange gaussian), a behaviour also found 
by T10. 

In Fig. [8] we show least-square fits to the [Z/H]-<j re- 
lationship in five narrow age-bins. These fits are indicated 
by the different colours with correspondingly coloured la- 
bels. The size of the bins were chosen to include a similar 
number of objects in each bin. Contours for two of the age 
bins are also shown, representing the fits with correspon- 
ding colours, and the grey-scaled filled contours represent 
the full sample. The fits to the age bins clearly show steeper 
slopes than for the entire red sequence population (dashed 
orange line). The fits to the age bins also show a parallel 
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Figure 7. Left hand panel: shows the relationship between the derived total metallicities and velocity dispersion. Rejuvenated objects 
with light-averaged ages smaller than 2.5 Gyr (see Fig. |6j are presented with cyan contours and the orange contours show the old red 
sequence population. The whole sample is shown as grey-scaled filled contours. The orange solid line is a least-square fit to the red 
sequence population and the parameters of the fit are given at the top of the panel. The dashed black line is the analogous fit from T10 
for comparison. Median lcr-errors are shown in the lower right corner. Right hand panel: shows the distribution of deviation in age from 
the least-square fit to the red sequence population with the same colour coding as in the left hand panel. The standard deviation of the 
gaussian fitted to the distribution of the red sequence population (orange line) is indicated in the upper right corner. 
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Figure 8. Contour plots for the relationship between total me- 
tallicity and velocity dispersion in bins of narrow log(age Gyr) 
intervals. Blue contours are for 0.6<log(age Gyr)<0.75 and red 
contours for 0.9<log(age Gyr)<1.0. Grey filled contours show the 
whole sample. Coloured lines are least-square fits to bins of vary- 
ing age intervals, indicated by the correspondingly coloured labels 
that also give the parameters of the fits. The orange dashed line 
is the least-square fit to the whole sample (see Fig. 



behaviour differing mainly in zero point offsets. This indica- 
tes a planar dependence on age and velocity dispersion for 
total metallicity, i.e. total metallicity correlate with velocity 
dispersion at fixed age and anti-correlate with age at fixed 
velocity dispersion. Su c h a be haviour has previously been 
found by iTrager etal] l|2000bl ). The large scatter in [Z/H] 
found for the entire sample is partly an effect of this planar 
dependence. The shallower slope for the entire red sequence 
sample is due to the generally older ages found for galaxies 
with a higher velocity dispersion and vice versa (see Fig. [6j . 



5.4 [Fe/H] 

Iron abundances derived using Eq. [T3]in Section [33] are pre- 
sented in the left hand panel of Fig.[9]as a function of velocity 
dispersion, with the same colour coding as in Fig. [U]for the 
results of this work. The least-square fit (orange line) shows 
no correlation between iron abundance and velocity disper- 
sion for the full red sequence population. Hence the [Fe/H]- 
cr relation is significantly flatter than the [Z/H]-<r relation- 
ship presented in th e previous sec t ion. Sim ilar patterns have 
also b een found by ITrager et al.l (|2000bl ) and I Graves et all 
(2007), but the slopes of their [Fe /Hj-cr relations ar e diffe- 
rent from zero. Similar to this work lPrice et al.l l|201ll ) find a 
relatively flat [Fe/H]-loga slope for a sample of 222 passive 
galaxies. 

The distribution around the fit to the red sequence po- 
pulation is shown in the right hand panel of Fig. [9] The reju- 
venated population (cyan gaussian) shows an offset (~0.15 
dex) towards higher [Fe/H] compared to the red sequence 
population (orange gaussian). This offset is more pronoun- 
ced than for [Z/H] (compare with the right hand panel of 
Fig. 01 • 

In Fig. [10] we show least-square fits to the [Fe/HJ-er 
relationship in five narrow age-bins, with the same colour 
coding and bin sizes as in Fig. [8] Similar to the total me- 
tallicity case, we find a planar dependence on velocity di- 
spersion and age, i.e. iron abundance correlate with velo- 
city dispersion at fixed age and anti-correlate with age at 
fixed velocit y dispersion. Such a behaviour ha ve again be- 
en fou nd by ITrager et alj (|2000bh and also bv lSmith et al.l 
(2009). The planar dependence is partly responsible for the 
large scatter in iron abundance found for the full sample. 
Compared to total metallicity we find flatter slopes for the 
[Fe/H]-er relationships in the age-bins as well as for the full 
sample. The fits to the age bins in Fig. [IJJclearly show stee- 
per slopes than for the full red sequence population sample 
(dashed orange line), since older ages are in general found 
for galaxies with a higher velocity dispersion and vice versa 
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Figure 9. Left hand panel: shows the relationship between iron abundance and velocity dispersion. Rejuvenated objects with light- 
averaged ages smaller than 2.5 Gyr (see Fig. [6j are presented with cyan contours and the orange contours show the old red sequence 
population. The whole sample is shown as grey-scaled filled contours. The orange solid line is a least-square fit to the red sequence 
population and the parameters of the fit are given at the top of the panel. Median lcr-errors are shown in the lower right corner. Right 
hand panel: shows the distribution of deviation in age from the least-square fit to the red sequence population with the same colour 
coding as in the left hand panel. The standard deviation of the gaussian fitted to the distribution of the red sequence population (orange 
line) is indicated in the upper right corner. 
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Figure 10. Contour plots for the relationship between iron ab- 
undance and velocity dispersion in bins of narrow log(age Gyr) 
intervals. Blue contours are for 0.6<log(age Gyr)<0.75 and red 
contours for 0.9<log(age Gyr)<1.0. Grey filled contours show the 
whole sample. Coloured lines are least-square fits to bins of vary- 
ing age intervals, indicated by the correspondingly coloured labels 
that also give the parameters of the fits. The orange dashed line 
is the least-square fit to the whole sample (see Fig. |9j 



5.5 [O/Fe] 

The relationship between [O/Fe] ratio and velocity disper- 
sion is presented in the upper left panel of Fig II II with the 
same colour coding as in Fig. [6] For the red sequence popu- 
lation we find a somewhat shallower slope for the [0/Fe]-<r 
relation (orange line) compared to the [«/Fe]-cr relation from 
T10 (dashed black line), but the general agreement is very 
good, i.e. a tight correlation between [O/Fe] and velocity 
dispersion. The shallower slope can be explained by metal- 
licity dependent C yields that depress the amount of O in 
more metal-rich systems (see discussion in Section \6. 31) . 

Delayed enrichment of Fe-like elements from SNIa ex- 
plosion compared to the prompt enrichment of cv-elements 
from SNII explosion, results in lower a/Fe ratios for objects 
with more extended star formation histories. Hence T10 in- 
terpreted the higher [cv/Fe] ratios found in more massive 
galaxies as shorter formation time-scales for such objects, 
and evidence for down-sizing of early-type galaxies. This in- 
terpretation gets support from the [0/Fe]-<r relation found 
in this work, since O belongs to the a-elements. 

The upper right panel of Fig. [TT] shows the distribution 
around the fit to the red sequence population. The reju- 
venated population (cyan gaussian) shows an offset towards 
lower [O /Fe] ratios compared to the red sequence population 
(orange gaussian). This offset is less pronounced compared 
to what T10 found for [a/Fe]. 



(see Fig. [6] and Section [5.3|) . Hence, at fixed age the effec- 
tive Fe yields are higher in more massive systems. Howe- 
ver, Fe enrichment is suppressed in older galaxies because 
of time-scale dependent contribution from SN la. Therefo- 
re the relation becomes flat for the overall population and 
the [Fe/H]-age anti-correlation is steeper than the [Z/H]-age 
anti-correlation. 



5.6 [Mg/Fe] 

The [Mg/Fe]-cr relationship is presented in the lower left 
panel of Fig. 1111 with the same colour coding as in Fig. [6] 
for the contours and least-square fit. For the red sequence 
population [Mg/Fe] clearl y increases with velocity dispers i- 
on, a pattern a lso fo und by Sanchcz-Bl azquez et al.l(l2006al). 
iGraves et all (|2007l ). ISmith et al.1 (120091 ) and iPrice etall 
(20T1|)- 

The [Mg/Fe]-<7 relationship shows a similar trend 
as [0/Fe]-cr (compare with the upper left panel of Fig. Ill[) . 
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Figure 11. Left hand panels: show the relationships for the derived O/Fe (upper) and Mg/Fe ratios (lower) with velocity dispersion. 
Rejuvenated objects with light-averaged ages smaller than 2.5 Gyr (see Fig. [6j are presented with cyan contours and the orange contours 
show the old red sequence population. The whole sample is shown as grey-scaled filled contours. The orange solid lines are least-square 
fits to the red sequence population and the parameters of the fits are given at the top of the panels. Median lcr-errors are shown in the 
lower right corners. Right hand panels: show the distribution of deviation in age from the least-square fits to the red sequence population 
for O/Fe (upper) and Mg/Fe (lower), with the same colour coding as in the left hand panels. The standard deviation of the gaussians 
fitted to the distributions of the red sequence population (orange lines) are indicated in the upper right corners. 



A steeper slope is found for the former relationship, which is 
in better agreement with the slope of [a/Fe]-<r relation from 
T10. Thus the T10 [a/Fe] ratios trace [Mg/Fe] rather than 
[O/Fc] as discussed in Section f5. II 

The distribution around the fit to the red sequence po- 
pulation is shown in the lower right panel of Fig. 1111 with 
the same colour coding as in Fig. [6] The rejuvenated popu- 
lation (cyan gaussian) is offset towards lower [Mg/Fe] ratios 
compared to the red sequence population (orange gaussian). 
This offset is more pronounced than for [O/Fe] (see upper 
right panel of Fig. 1111 Section [5.5[) . This is again in better 
agreement with the results for [a/Fe] in T10, which sug- 
gests that the [a/Fe] ratio derived in T10 is closest to the 
true [Mg/Fe] ratio. 

5.7 [C/Fe] 

The [C/Fe]-cr relationship is shown in the upper left panel of 
Fig. 1121 with the same colour coding as in Fig.[6]for the con- 
tours and the least-square fit to the red sequence population. 
The [Mg/Fe] -<j relationship is included for comparison (das- 
hed orange line). Considering the red sequence population 



we find a very similar trend for the [C/Fe]-er relationship (so- 
lid orange line) as for [Mg/Fe]-<r (dashed orange line), but 
the latter trend is slightly steeper. The stronger trend for 
[C/Fe] can be explained by metallicity dependent C yields 
(see discussion in Section f6 . 3 . 

Increasing [C/Fe] ratios with velocity dispersi- 
on has also been found by ISanchez-Blazquez et al.l 

d |Sn 



1 2006 j). [Graves et al l (120071) and ISmith et all ( 2009 



ISanchez-Blazquez et al.l l|2006al ) report a stronger increase 
in [Mg/ Fe] than [C/Fe] fo r a sample of 98 ear ly-type 
galaxies. iGraves et ail (|2007l ) and lSmith et al.l l|2009h study 
~6000 red sequence SDSS early-type galaxies and 147 red 
sequence galaxies, respectively, and find, in agreement with 
this work, a shallower trend for [ Mg/Fe] than for [C /Fe]. A 
similar pattern was also found bv lPrice et alj (|201lf ). 

The distribution around the fit to the red sequence po- 
pulation is shown in the upper right panel of Fig. 1121 The 
rejuvenated population (cyan gaussian) shows an offset of 
lower [C/Fe] ratios compared to the red sequence populati- 
on (orange gaussian). Again very similar to [Mg/Fe] we find 
a clear offset towards lower [C/Fe] ratios for the rejuvenated 
population compared to the red sequence population. Hence 
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Figure 12. Left hand panels: show the relationships for the derived C/Fe (upper) and N/Fe ratios (lower) with velocity dispersion. 
Rejuvenated objects with light-averaged ages smaller than 2.5 Gyr (see Fig.0 are presented with cyan contours and the orange contours 
show the old red sequence population. The whole sample is shown as grey-scaled filled contours. The orange solid lines are least-square 
fits to the red sequence population and the parameters of the fits are given at the top of the panels. Median lcr-errors are shown in the 
lower right corners. Right hand panels: show the distribution of deviation in age from the least-square fits to the red sequence population 
for C/Fe (upper) and N/Fe (lower), with the same colour coding as in the left hand panels. The standard deviation of the gaussians 
fitted to the distributions of the red sequence population (orange lines) are indicated in the upper right corners. 



we find that the [a/Fe] ratio derived in T10 resembles both 
[Mg/Fe] and [C/Fe]. In fact the slope of the [a/Fe]-tr re- 
lationship from T10 (0.33) falls right in between the slope 
for [Mg/FeJ-a (0.30) and [C/Fe]-ff (0.35). Following the dis- 
cussion in Section 15.11 this is probably due to the fact that 
T10 use absorption line indices sensitive to various element 
abundance variations to derive [a/Fe] without considering 
individual abundance ratios. These indices show the stron- 
gest signals to Mg and C, such that their derived [a/Fe] 
ratios reflect [Mg/Fe] and [C/Fe] rather than [O/Fe]. 

In Section 16.31 we discuss the controversy of the origin 
of C and that the similarity between the derived [Mg/Fe] 
and [C/Fe] ratios reflect a significant contribution of C from 
massive stars. 

5.8 [N/Fe] 

The [N/Fe]-cr relationship is shown in the lower left panel 
of Fig. 1121 with the same colour coding as in Fig. [6] for the 
contours and least-square fit to the red sequence population. 
The [Mg/Fe]-i7 relationship is shown for comparison (dashed 
orange line). Despite a fairly large scatter we find a clear cor- 



relation between [N/Fe] and velocity dispersion for the red 
sequence population, i.e. [N/Fe] strongly increases with ve- 
locity dispersion. This correlation shows a significantly stee- 
per trend than the analogous found for [Mg/Fe] and [C/Fe]. 
The [N/Fe] ratios are also offset towards lower abundance 
ratios by ~0.2 dex compared to [Mg/Fe] and [C/Fe]. An in- 
cr ease in [N/Fe] with velocity disp e rsion has also been found 
bvlSanchez-Blazquez et all |2006al ). lGraves et al.l J2007f) and 
ISmith et alj ""(f2009). The result of ISanchez-Blazquez et al.l 
(2006a) show steeper trends w ith velocity disper sion for 
[Mg/Fe] than for [N/Fe], while ISmith et all |2009l ) report 
similar trend s, steep er th an [Mg/Fe], f or [N/ Fe] and [C/Fe]. 
iGraves et al. | (|2007h and iPrice et al.l (|201ll ) find, again in 
agreement with this work, the steepest trends with veloci- 
ty dispersion for [N/Fe]. Also in agreement with this work 
they find [N/Fe] to be offset towards lower abundance ra- 
tios compared to [Mg/Fe] and [C/Fe]. However, this offset 
is dependent on their adopted, fixed value for [O/Fe] (see 
Section |4~4)) . 

The distributions around the fit to the red sequence 
population are shown in the lower right panel of Fig. [T2l The 
rejuvenated population (cyan gaussian) shows no clear offset 
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Figure 13. Left hand panels: show the relationships for the derived Ca/Fe (upper) and Ti/Fe ratios (lower) with velocity dispersion. 
Rejuvenated objects with light-averaged ages smaller than 2.5 Gyr (see Fig. [6j are presented with cyan contours and the orange contours 
show the old red sequence population. The whole sample is shown as grey-scaled filled contours. The orange solid lines are least-square 
fits to the red sequence population and the parameters of the fits are given at the top of the panels. Median lcr-errors are shown in the 
lower right corners. Right hand panels: show the distribution of deviation in age from the least-square fits to the red sequence population 
for Ca/Fe (upper) and Ti/Fe (lower), with the same colour coding as in the left hand panels. The standard deviation of the gaussians 
fitted to the distributions of the red sequence population (orange lines) are indicated in the upper right corners. 



compared to the red sequence population (orange gaussian), 
considering the fairly large scatter in the data. 

The differences in the derived [C/Fe] and [N/Fe] ratios 
are interesting as both elements are expected to be enriched 
through the same processes. A thorough discussion on the 
origins of C and N are given in Section T6.3I and [6.41 respec- 
tively, and related to the derived element ratios. 



5.9 [Ca/Fe] and [Ti/Fe] 

Trends with velocity for [Ca/Fe] and [Ti/Fe] are shown in the 
upper left hand and lower left hand panels of Fig. [l3l respec- 
tively, with the same colour codings as in Fig. [6] The derived 
[Ca/Fe] ratios show very weak trends with velocity dispersi- 
on and a fairly small scatter. The trend is significantly flatter 
than compared to the [Mg/Fc]-ir relationship (dashed orange 
line) and [Ca/Fe] is generally lower than [Mg/Fe] by ~0.1- 
0.2 dex. This implies an underabundance of Ca compared 
to Mg, such that Ca scales more closely with Fe rather than 
Mg. The [Ti/Fe] ratios show a very large scatter, which is 
due to just a weak response to Ti found for the Lick indices 
(see Fig. [TJ, requiring very high S/N observations to cons- 



train [Ti/Fe]. Thus, we can't draw any strong conclusions on 
the [Ti/Fe]-cr derived trends. Still, when compared to the 
[Mg/Fe]-cr relationship (dashed orange line) the data sug- 
gest an overall under-ab undance of Ti compared to Mg with 
an offset similar to Ca. iMilone. Barbuv k, Schiavonl (|2000h 
instead find that [Ti/Fe] « [Mg/Fe] using the TiO bands for 
a sample of 12 galaxies. 

The under-abundance of Ca compared to Mg, C 
and O confi r ms p r evious findings (| Thomas et al. 2003b; 



and U confi r ms p r evious h ndmgs ( 1 nomas ct al. 20031 
iGraves et all 120071 : ISmith et al.1 120091 : iPrice et all 12011 



This has been interpreted as a contr ibution of Ca from SNIa 
besides SNII (|Thomas et al.1 1201 lal ). Since Ti is a heavier 
element than Ca the under-abundance of both of these ele- 
ments imply that type la Supernovae also contributes to Ti. 
These results hint that the contribution from type la Super- 
novae is dependent on atomic number for the a-elements. 
This pattern of lower Ca and Ti abundances compared to 
Mg, O, and C has also been found for the s tellar populati- 
ons of the Milky Way (|Thomas et al.lr2011al . and references 
therein). Hence heavy a-elements are universally produced 
in SNIa. 

The right hand panels of Fig. [13] show the distributions 
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Table 4. Average offsets in abundance ratios, metallicity and 
iron abundance for the rejuvenated population compared to the 
red sequence population. First column states the parameters, se- 
cond column the offsets and the third column gives the standard 
deviation of the distributions for the rejuvenated population. 



Parameter 


offset 


a 


[O/Fe] 


-0.0265 


0.129 


[C/Fe] 


-0.0577 


0.0931 


[N/Fc] 


-0.0526 


0.183 


[Mg/Fe] 


-0.0539 


0.0980 


[Ca/Fe] 


-0.0600 


0.130 


[Ti/Fe] 


-0.143 


0.360 


[Z/H] 


0.106 


0.248 


[Fe/H] 


0.147 


0.266 



around the red sequence populations for Ca (upper panel) 
and Ti (lower panel), with the same colour coding as in 
Fig. [U The rejuvenated population (cyan gaussian) show 
weaker abundance ratios compared to the red sequence po- 
pulation. This is however less pronounced than for Mg and 
Ca because of a significant contribution from both SN II 
and SN la to Ca. Clearly the standard deviation of the reju- 
venated population is fairly large indicating that the exact 
position of the rejuvenated [Ca/Fe] peak is not very well 
defined. For Ti the peak of the rejuvenated population is 
even less well defined and no conclusions regarding the this 
population can be drawn from Ti. 



6 DISCUSSION 

We present a method for simultaneously deriving the ele- 
ment abundance ratios [C/Fe], [N/Fe], [Mg/Fe], [Ca/Fe] 
and [Ti/Fe] for unresolved stellar populations, together with 
the classical stellar population parameters age, [Z/H] and 
[O/Fe] (representing [a/Fe], see Section [33}. The method 
is based on the new flux-calibrated stellar p opulation mo- 
dels of absorption line indices presented in iThomas et al.l 
(|2011bl). It has been cali brated on galactic globular clusters 



Thomas et al.l (|2011al ). showing a good agreement for the 



derived parameters with the corresponding of high resoluti- 
on spectra of individual cluster stars. We apply the method 
to a sample of 3802 SDSS early-type galaxies. The same 
sample of galaxies was used in T10 who derived the stellar 
populations parameters age, [Z/H] and [a/Fe], based on the 
TMB/K models. Hence the analysis presented here can be 
regarded as an extension to the work of T10. The results 
derived here are in good agreement with the results of T10, 
i.e. age, total metallicity and [a/Fe] increase with increa- 
sing velocity dispersion. This is remarkable given that we 
use new models and a new method for measuring the stel- 
lar population parameters that includes individual element 
abundance ratios. 



6.1 Rejuvenated population 

T10 find the sample of visually identified early-type gala- 
xies to have, in addition to a dominant red sequence gala- 
xy population, a sub-population of rejuvenated galaxies, i.e. 



galaxies with mainly old stellar populations that have ex- 
perienced minor recent star formation producing stars that 
overshine the dominant old stellar populations and thus mi- 
micking overall young stellar populations. The rejuvenated 
galaxy population was identified by a secondary peak in the 
derived age distribution showing younger ages than the pri- 
mary peak and consists mainly of low mass early-type ga- 
laxies. This was supported by higher total metallicities and 
most importantly by lower [a/Fe] ratios and detection of 
residual star formation through presence of emission lines. 
Schawinski et al. (2007) and T10 present the emission li- 
ne diagnostics and show that most of the galaxies in the 
the young subpopulation have emission lines caused by star 
formation activity. 

The higher total metallicities imply that residual star 
formation in the rejuvenated galaxies occurred in an ISM 
containing a fraction of gas that have been chemically en- 
riche d compared to pristine metal-poor gas (Th omas et alj 
l2010l y In this work we further show that the rejuvenated 
population shows a stronger excess in iron abundance com- 
pared to total metallicity. This is due to the delayed Fe en- 
richment from SN la, which consequently produces lower 
[a/Fe] ratios. Hence the extended star formation histories 
of the rejuvenated galaxies allow enough time for the ISM 
to be highly enriched in Fe. 

Table 2] summarises the offsets between the rejuvenated 
and red sequence population derived in Section 15.21 - 15.91 
The standard deviations (a) for the distributions of the re- 
juvenated population are also included. The result of T10 
is reproduced in this work, where we find younger ages to- 
gether with higher [Z/H], higher [Fe/H] and lower [E/Fc] 
ratios for the rejuvenated population. In more detail we find 
[O/Fe], [Mg/Fe], [C/Fe], [N/Fe] and [Ca/Fe] to be offset to- 
wards lower abundance ratios by 0.03-0.06 dex. However, 
the precision of the offsets varies significantly and shows the 
highest accuracy for [Mg/Fe] and [C/Fe], while the offset for 
[N/Fe] is not well determined. [Ti/Fe] shows the largest off- 
set, but to a very low precision. The accuracy of the offsets 
is due to the sensitivity of the indices to the variation of the 
different element abundances (see Section |3.4[1 . [Z/H] and 
[Fe/H] show positive offsets for the rejuvenated population, 
with [Fe/H] being more offset by ~0.04 dex, i.e. reflecting 
the offset found for [E/Fe]. 

6.2 Environment 

T10 find environmental dependencies for the fraction of reju- 
venated galaxies, while the red sequence galaxy population 
is unaffected by the environment for the classic parameters 
age, total metallicity and a/Fe. In brief the fraction of reju- 
venated galaxies increases with decreasing velocity disper- 
sion and environmental density. We reproduce these results 
and do not further discuss environmental dependencies for 
the classic parameters, the reader is instead referred to T10 
for more details. We have also looked at environmental de- 
pendencies for the various element abundance ratios and no 
such can be seen. 

Hence we do not co nfir m earlier findings of 
ISanchez-Blazquez et al.l l|2003f ) and ISanchez-Blazquez et al.l 
l|2006al ). who report over-abundances of carbon and ni- 
trogen in environmen ts of lower density fo r 98 early-type 
galaxies. The work of lClemens et al. (2006) studies a large 
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Figure 14. [C/O] (black contours) and [C/Mg] (orange contours) 
as functions of velocity dispersion for the red sequence populati- 
on. The solid lines are least-square fits to the relationships with 
corresponding colours and the dashed black line indicate solar 
abundance ratios. The slope and intercept of the least-square fits 
are given by the labels. 

sample of SDSS (Sloan Digital Sky Survey. I York et al]l200d ) 
early-type galaxies and find, in agreement with the results 
derived here, that the environment does not influence the 
enhancement of carbon. 

6.3 Lower limit on formation time-scales from 
[C/Mg] 

Fig. ll4l shows contours and least-square fits for [C/O] (black 
contours and lines) and [C/Mg] (orange contours and lines) 
as functions of velocity dispersion for the red sequence po- 
pulation. Both of these abundance ratios are close to solar 
values. However, differences are noticeable. [C/Mg] is ~0.0 
for the most massive systems and ~-0.05 for the less mas- 
sive galaxies. [C/O] is instead ~0.05 for the most massive 
systems and ~-0.05 for the less massive galaxies. The sources 
of C production have bee n the subject of a long standing de- 
bate. iHenrv et al] (2000) compared abundances of galactic 
and extra-galactic HII regions to models of chemical evo- 
lution. To match the data they needed stell ar yields where 
massiv e stars dominate the production of C. ICescutti et al.l 
(2009) found that stellar yields with metallicity dependent 
C produced in massive stars are needed to match models of 
chemical evolution with element abundance ratios derived 
from individual stellar spectra of the galactic stellar popula- 
tions. Stellar yields take mass-loss and rotation into account. 
The latter cause newly synthesised C to be brought to the 
surfa ce layers and ejected into the ISM l|Mevnet fc Maederl 
2002) by mass-loss. The mass-loss rate is dependent on me- 
tallicity such that more metal-rich massive stars experience 
higher mass-loss rates and consequently contribute higher 
abundances of C. 

Hence recent results indicate that massive metal-rich 
stars are a major C source. Still, intermed iate mass stars 
do al s o contribute significant amounts of C (jRenzini fc Volil 
Il98ll ; I van den Hoek fc Groenewegenl 1 19971 ). Dredge- up on 
the asymptotic giant branch (AGB) phase of stellar evolu- 
tion bring C up to the surface layers. This occurs following 
He-shell flashes when the convective envelope reach down 



to the inner layers where He-burning has taken place. The 
dredged-up C is then ejected i nto t he ISM by mass-loss. 
Ivan den Hoek fc Groenewegenl l|l997l ) show that C yields 
from 3Mq stars are significant. To reach solar C/Mg values 
star formation must continue over long enough time-scales 
to allow for the contribution of C from both massive and 
intermediate mass stars. This sets a lower limit for the star 
formation t ime-scales of ~0-4 Gyr, which is the lifetim e of a 
3M stars (|Castellani et al.lll992l : iBertelli et al.ll2009T ). This 
is thus the lower limit for the formation of the most massive 
early- type galaxies of our sample, since these have [C/Mg]~0 
(see Fig. ll4p . In a scenario where self-enriched bursts of star- 
formation build up the stellar populations of early-type ga- 
laxies, the bursts must last for at least ~0.4 Gyr. However, 
the formation time-scales must be long enough to allow the 
enrichment of the ISM to produce the high metallicities ob- 
served in massive early-type galaxies. Indeed, a similar time- 
scale (~0.4 Gyr) for the onset of galactic winds in models 
of chemical evolution can produce massive early-type ga- 
laxies with metallic ities in agreement with observed values 
l|Pipino et al.ll2010l ). 

Although weak, the trend of increasing [C/Mg] with 
velocity dispersion (see Fig. [T4|l could be caused by a metal- 
licity dependent production of C in massive stars. The more 
massive galaxies are also more metal-rich (see Section [5.31) . 
Maeder (1992) predict metallicity dependent C yields from 
increasing mass-loss in more metal-rich, massive stars. The 
net result is an ejection of large amount of Carbon before 
this element is turned into heavier elements. The side ef- 
fect is that less Carbon is available for producing Oxygen. 
Thus Maeder (1992) also predicts a depression of Oxygen 
along with the enhancement of Carbon at higher metalli- 
cities, while at low metallicities the opposite situation is 
apparent. This behaviour was confirmed by e.g. Cescutti 
et al. (2009) for the Milky Way. In Fig. [H] we have also 
seen that [C/O] increases with velocity dispersion such that 
super-solar and sub-solar [C/O] ratios are found for the most 
massive and least massive galaxies, respectively. Hence the 
flatter slope found for the [O/FeJ-tr relationship compared 
to [Mg/Fe]-<7 can be explained by the metallicity dependent 
balance between C and O yields from massive stars, since 
higher velocity dispersion galaxies are more metal-rich. The 
effects discussed above are expected to be weak considering 
the short metallicity range covered. 

6.4 The Nitrogen puzzle 

The origin of N is probably even more debated than the ori- 
gin of C. It remains controversial whether N is a secondary 
and/or a primary element. N is produced in the CNO-cycle 
from C and O. Primary N then comes from the CNO pro- 
cessing in stars, while secondary N is produced from C and 
O already present in the c ollapsing gas-clouds forming the 
stars (e.g. Matteucci 198jJ). The amount of primary N pro- 
duced is proportional to the abundances of other primary 
elements (e.g. C and O). The N abundance from secondary 
production is instead proportional to the initial C and O 
abundances. Primary and secondary production in interme- 
diate mass stars, with masses in the range 4M0<M<8M0, 
are believed to be the dominant sources of N, where dredge- 
up, hot bottom burning and mass-loss eject N into the ISM 
|Renzini fc Volilll98ll ; Ivan den Hoek fc Groenewegenlll997r i . 
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iMevnet fc Maederl (|2002h show that N can also be ejected by 
massive stars through the inclusion of stellar rotation. Ob- 
servations indicate that primary N in low metallicity massi- 
ve star s is required to e x plain observed trends of abundance 
ratios. Ilzotov fc Thuanl (|l999h found HII regions in low me- 
tallicity blue compact galaxies to show a flat log(N/0) trend 
with 12+log(0/H) for their lowest metallicity galaxies, while 
the log(N/0) ratio starts to increase at a certain metallicity 
(7.6<12+log(0/H)<8.2). The interpretation is that the pri- 
mary N is produced by massive stars at low metallicities and 
the contribution from intermediate mass stars is delayed in 
time and kicks in at higher metallicities. 

Fig.lT5lshows the [N/0]-er (black contours) and [N/Mg]- 
a relationships (orange contours) together with least-square 
fits (solid lines) with corresponding colours. Fit parameters 
are given by the labels. An under-abundance of N of ~0.2 
dex compared to O and Mg is accompanied by a significant 
slope that is more prominent for the [N/0]-<r relation. The 
different origins of N compared to Mg and O make the N/Mg 
and N/O ratios useful formation time-scale indicators. Fol- 
lowing our discussion of C/Mg as a lower limit time-scale 
indicator (see Section [57T)) . the low N/Mg ratios suggest that 
the formation time-scales of the red sequence sample are too 
short for the full N production to be reached. This could be 
the case if low mass stars with long life-times contribute 
significantly to t he production of N. Such indications have 
been reported bv lThuanT Pilvugin fc Zinchenkd l|2010D from 
emission lines of SDSS star-forming galaxies, who find that 
significant amounts of N are produced in 1.5-2Mq stars with 
lifetimes of 2-3 Gyr. 

In the low velocity dispersion regime (of our sample) we 
find [N/Fe]<0.0 (see lower left panel of Fig. [12]). This would 
require N to be produced over longer time-scales than Fe 
and could constrain the upper formation time-scale limit for 
the low mass systems, if the production sites of N were well 
co nstrained. Considering the ste l lar m ass range suggested 
by Eh uan, Pilvugin fc Zinchenkd (|2010h for N contribution, 
an upper limit of ~2 . 5 Gyr (life-time of a 1.5 Mq stars, 
ICastellani et al.lll992l ; iBertelli et al.l I2008T ) may be implied 
for the galaxies in our sample with 1.8< a <2.0 that have 
[N/Fe]<0.0. Future models with better constrained N yields 
are needed to quantify this to a higher degree. 

However, besides the under-abundance of N the slopes 
of the N/Fe-er, N/Mg-cr and N/O-ct relationships suggest 
that abundance ratios are higher in more massive galaxies. 
Following the discussion above longer formation time-scales 
would be required for such systems compared to the lower 
mass galaxies. This is obviously in contradiction to the shor- 
ter formation time-scales of more massive systems implied 
by the better constrained, in terms of the sources of Mg 
production, Mg/Fe ratios. Hence in the scenario discussed 
above the N/E ratio does not hold as formation time-scale 
indicator for the most massive systems 

Instead the Mg/Fe ratio would set an upper limit on the 
formation time-scales of the most massive systems and N/Fe 
provides an additional upper limit on lower mass galaxies. 
On the other hand if N is mainly produced in stars with 
masses above 4Mq , the delayed enrichment of Fe from SNIa 
will result in a correlation with velocity dispersion also for 
[N/Fe] (see discussion in Section ft>.3[) . Still this correlation 
is stronger than for [O/Fe], [Mg/Fe] and [C/Fe]. In either 
case, the higher N/E ratios in the more massive systems 



4 [ 

slope/intercept: 
■ [N/0] 0.89±0.02/-0.77±0.0S 

0.2 " [N/Mg] 0.17±0.02/-0.S5±0.04 
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Figure 15. [N/O] (black contours) and [N/Mg] (orange contours) 
as functions of velocity dispersion for the red sequence populati- 
on. The solid lines are least-square fits to the relationships with 
corresponding colours and the dashed black line indicate solar 
abundance ratios. The slope and intercept of the least-square fits 
are given by the labels. 



still need to be explained. Independent of the lower stellar 
mass limit for significant N contribution, at least three diffe- 
rent scenarios can steepen the N/Fe-cr, N/O-er and N/Mg-a 
trends 

i) Metallicity dependent secondary N producti- 
on. As mentioned secondary production of N will cause a 
strong increase in N abundances with increasing metallicity. 
With the higher velocity dispersion galaxies being overall 
more metal-rich the steep [N/FeJ-er trends could intuitively 
be explained by the secondary production of N in interme- 
diate mass stars. The metallicity range covered is relatively 
small, however, and it seems contrived if metallicity effects 
in secondary N production can explain the steep [N/Fe]-<r 
relationship alone. 

ii) Primordial gas inflow. The ratios between secon- 
dary and primary elements are useful indicators of primor- 
dial gas accretion, while the ratios bet ween two primary ele- 
ments are unaffected by such events (jKoppen fc Edmunds! 
1 19991 ; iHenrv et al.|[2000l ). Hence N/ Mg and N/O ratio s are 
tracers of primordial gas inflow (e.g. lHenrv et al" 2000) that 
results in lower N abundances compared to Mg and O. 
Fig. ll5l shows the [N/OJ-ct relationship. An under-abundance 
of N of ~0.2 dex is accompanied by a significant slope. This 
trend could be caused by a stronger dilution of the ISM 
in low mass galaxies. If the high velocity dispersion gala- 
xies experienced very intensive gas accretion over time-scales 
shorter than the star formation time-scale, high N abun- 
dances are achieved since dilution of the ISM does not take 
place over the entire star formation history. If low velocity 
dispersion galaxies instead experience gas accretion along 
with star formation the ISM instead gets diluted with pri- 
mordial gas over the entire star formation history and the 
N abun dances become lowe r . This scenario fits with the re- 
sults of iDekel fc Birnboim I (2006) who show that the shut- 
down of gas supply from cold streams is dependent on halo 
mass, such that it continues over longer time-scales in less 
massive halos. Such a scenario can steepen the slope of the 
[N/0,Mg]-<r relationships and consequently [N/Fe]-<r, but al- 
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so for this case it is not clear if it can produce the observed 
trends alone. 

iii) Varying initial mass function (IMF). IMF 

variations have been considered throughout the litera- 
ture to explain variations of element ratios in early- 
type galaxies, especially the tr end between [ Mg/Fe] 



and velocit y 
Trager et all feoOOb : 



dispersion 



(e.g. IThomas et all [l999bl; 



Maraston et all [20031: iMatteuccil 
1994 ISanchez-Blazquez et alj 12006 j ; ISmith et all I2009T ). 

However, the variations considered a l so affe ct absorption 
line indices dSanchez-Blazquez et al. 2006a ) and galaxy 
scalin g relations ( Renzini fc CiottJ " 19931 ; Maraston et al.l 
120031 . and references therein) in negative ways questio- 
ning the plausibility of IMF dependencies on galaxy 
mass. Since intermediate mass stars are believed to be 
the dominant sources of N, IMF variations with galaxy 
mass dependent weights in the (approximate) stellar 
mass range 5Mp;<M<8Mp) (partic ularly high N yields, 
Ivan den Hoek fc Groenewegenl 1 19971 ) are needed to vary 
the N abundances. IMF variations with emphasis on such 
a specific mass range have not been observed and are 
unlikely to occur. This would affect C abundances as well, 
hence the more or less constant C/Mg ratios (over the 
velocity dispersion range covered) does not favour the IMF 
variations discussed. 

The first two mechanisms (i-ii) could work together 
to cause the observed trends with velocity dispersion. To 
evaluate if this is plausible simulations of chemical evo- 
lution are n eeded taking these mechanisms into account. 
However, in IPipino et al.1 |2010l ) (P10) up-to-date models 
of chemical evolution fail at reproducing observed [N/Fe]- 
mass trends of early - type g alaxies. P10 adopt the models of 
IPipino fc Matteuccil (2004 ) (PM04) that reproduce the ob- 
served pattern of increasing Mg/Fe ratios with increasing 
stellar mass in early-type galaxies, by assuming star forma- 
tion histories compatible with the down-sizing scenario and 
Mg contributed by type II Supernovae. The models also re- 
produce the observed abundances of Ca by adopting yields 
where contributions to this element come from both type II 
and type la Supernovae. 

Implementing different recipes of stellar yields P10 com- 
pare modified versions of the PM04 model s to th e observed 
element abundance ratios of [Graves et al.1 (|2007I i. The mo- 
dified versions mainly differ in the prescription of stellar 
mass-loss and rotation, which impacts on the abundances 
of C and N. Stellar rotation causes a mixing of elements in 
different stellar layers and mass-loss the ejection of newly 
synthesised element into the ISM. With the new prescripti- 
ons P10 are able to match the observed [C/Fe]-<r trend, both 
in zero-point and slope. The results indicate that there must 
be a substantial contribution to the production of C from 
massive stars and metallicity dependent yields possible due 
to metallicity dependent mass- loss rates. For [N/Fe] they in- 
stead find a very large scatter between the different models 
such that the observed steep trend of the [N/Fe]-cr relati- 
onship can not be simultaneously matched with the overall 
high [N/Fe] ratios. Also, to reach these high abundance ra- 
tios they must adopt a prescription that is not physically 
justified. 



7 CONCLUSIONS 

We present light-averaged ages, metallicities and element 
abundance ratios for 3802 SD SS early-type galaxies drawn 
from the MOSES catalogue l|Schawinski et al.1 [20071 ) with 
visual morphology classifications. Using the flux-calibrated 
TMJ models of absorption line indices, which are based 
on the MILES stellar library, we have developed a me- 
thod for simultaneously deriving the element abundance ra- 
tios [C/Fe], [O/Fe] (inferred from [a/Fe]), [N/Fe], [Mg/Fe], 
[Ca/Fe] and [Ti/Fc]. The models are well calibrated with 
galactic globular clusters with independent measurements 
of stellar population parameters and element ratios (TMJ, 
IThomas et"al]|2011al ). 

We study the relationships between the stellar popu- 
lation parameters and galaxy stellar velocity dispersion. In 
agreement with the literature stellar population age and to- 
tal metallicity correlate with velocity dispersion. [Fe/H] in- 
stead does not show such a correlation over the entire pa- 
rameter range covered, but for a fixed age a steep trend is 
found for the [Fe/H]-<r relation. This trend is shallower than 
the analogous for [Z/H] due to suppressed Fe enrichment 
in more massive galaxies because of time-scale dependent 
contribution from SN la 

Similar trends are found for [O/Fe], [Mg/Fe] and [C/Fe], 
i.e. strong correlations with velocity dispersion in agreement 
with the literature. The first two are expected to be similar, 
since both O and Mg belong to the group of a-elements 
produced in massive stars through type II Supernovae. This 
is also in favour of the down-sizing scenario of early-type 
galaxies that set an upper limit on the star formation time- 
scales and wher e more massive syst ems experience shorter 
time-scales (e.g. IThomas et aT[|2010h . 

The C/Mg ratios are close to solar values, which instead 
sets a lower limit for the formation time-scales of early-type 
galaxies. Stars with masses down to ~3 Mq contribute signi- 
ficantly to the production of C. To reach solar C/Mg ratios 
formation time-scales need to be long enough for such stars 
to eject C into the ISM. The inferred lower formation time- 
scale limit is then ~0.4 Gyr, which is the life-time of a 3 Mq 
star. 

The [N/Fe] ratios are overall lower by ~0.2 dex compa- 
red to [O/Fe] and [Mg/Fe] and the trend with velocity di- 
spersion is very steep, i.e. more massive galaxies have signi- 
ficantly higher [N/Fe] ratios. The observed [N/Fe]-<r trends 
are difficult to interpret due to uncertainties in the origin of 
N. The zero-point and slope of this relationship can not be 
simultane ously matched by up-to-date models of chemical 
evolution (|Pipino et al.ll20ld ). Either the theoretical stellar 
yields have to be increased by a significant factor or other 
prescriptions have to be incorporated into the models that 
affect the N yields. Such prescriptions could be: 1. N yields 
with a stronger dependence on metallicity, since more massi- 
ve early-type galaxies are more metal-rich. 2. A dependence 
on galaxy mass for the ratio between the time-scale of star 
formation and the time-scale of primordial gas inflow, which 
affects the N/O, N/Mg and N/Fe ratios due to the secondary 
nature of N. 

We do not find any dependence on environmental den- 
sity for the element ratios studied. This is in contradiction 
to previous studies that have reported environmental depen- 
dencies for C and N abundances. Hence difference formation 
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scenarios for field and cluster early-type galaxies can not be 
inferred from the element ratios studied in this work. 

The [Ca/Fe] ratios do not correlate significantly with 
velocity dispersion and are close to solar values over the 
entire velocity dispersion range covered. Although tentati- 
ve, due to large errors, Ti shows a behaviour similar to Ca. 
This indicates an atomic number dependent contribution 
from type la Supernovae to the production of Q-elements, 
i.e. the yields from type la Supernovae are higher for heavier 
a-elements. This is now universally found since similar pat- 
terns have been found in th e stellar populations of the Milky 
Way (|Thomas et alJlioilal . and references therein) and puts 
strong constraints on supernova nucleosynthesis. 
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